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CHAPTER^ 1 

The upsuf^rge of our indutrial civilization in the 
first half of the twentieth century was founded upon 
fossilfuels- coal and oil. But there are many "have-not" 
nation^ and even contries rich in this fuels are now seeing 
rapid depletion of their reserves At this juncture uranium 
has come to the rescue as a hope for the future. The worlds 
uranium and thorium, it is estimated, represent an energy 

l£. 

reservs/somewhere between ten and hundred times larger than 
coal. Even so fissinable fuels t®©- are an exhaustible 
supply. At the rate at which the world's energy needs are 
expanding practically all the econcmically recoverable 
uranium could as coal, might be exhausted within another 
century 'or so. 

Besides limited uranium reserves, fission power 

also presents a more immediate problem; namely, disposal 
& 

of its radiactive wastes. 

All this helps to explain the drive, if not race, 
to f^wid out whether thermonuclear power can be tapped and 
put to work. If the fusion feaction can be made to yield 
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power f it will solve forever both the fuel supply 
problem and the problem of radioactivewastes . The basic 
fusion element deuterium, is as inexhaustible as the 
oc€ans, and fusion produces no appreciable amount of 
radioactive by products. 

Nuclear fusion is not exactly a new phenomenon . 

It has been generating the power of the Sun and other stars 
for billions of years. But to create and control fusion power 
on earth is a problem of a totally different order from 
harnessing fission. It is undoubtedly the most difficult 
project ever presented to scientists and engineers. 

The fission releases energy, because part of the 

mass of the fusing nuclec ?s transformed into energy 

2 . , 

according to Einstein's famous equation E= me . ‘A 
continuous fusion reaction, is analogous to the familiar 
process of combustion. To stick together or fuse, the molecules 
must collide violently, which means the material must be 
heated. Thrree condition are needed to burn a chemical fuel 
and harness its heat to do work; (1) The fuel must be 
rais^ to its ignition point. (2) There must be enough of 
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fusion reactions have sioidiif icant reaction cross-=section 
at relative energies ranging betv^een ten keV and hundred 
keV, At these temperature the matter is completely 
ionized and exists in the form of plasma in which the 
ion and electron densities are equal to maintain charge 
neutrality. 

There are basically two approaches to fusion. The 
first of these is magnetic confinement of plasma. The reacting 
plasma is too hot to be contained by materialwallS/ and it is 
therefore proposed to use a "magnetic bottle" i.e a magnetic 
field configuration which curls up the particle orbits in 
such way as to greatly reduce their escape rate from the 
reaction volume. Speaking in a very gen oral way y the confined 
state of plasma is far from equilibrium, and there are t 
many mechanismsleading to the break up, or dispersal, of the 
confined plasma. Research over many years he s lead to an 
understanding of these escape rautes and many have been 
eliminated. There is , however one phenomenon which is 
not yet understood; it was discovered by Bohm in 1949, 
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According to him, there exists an amomalous diffusion 
process in a magnetic field B. The particles diffuse 
with a velocity given by 

v„= - ckT cr V, 8.64xl0^T 

° "e - "e 

where Yg is a factor which theory cannot account for in a 
straightforward way and which according to Bohm is 16. 

Much research in recent years has consisted of attempts 
to discover XA/hich circumstances are reponsible for this 
anomalous diffusions, presumably caused by turbulent 
fluctuation process. 

Putting, = n^r 

shows that Bohm diffusion limits the life time, T = r/Vj^, 
of a plasma confined in cylindrical geometry with a radius 
r , to a value approximately given be 

T = 10^ ("Yg B)/T 

where"/, is an improvement factor which may be raised by 

X A 

controlling the fluctuation level. There B is measured in 

8 

Gauss and T in Kelvin so that for T = 10 and a field of 
hundred kG and r = 100 cm, T= 16ms , if = 1 . 
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Magnetic confinement systems are characterized 
by the ratio of thermal energy 2nkT and magnetic energy 

bVs , 

^ = 16 A" n kT/B^ 

The parameters of a possible reactor can be 

determined if we also consider the power loading of the 

reactor 'Halls. A figure of 500W/cm is commonly taken as a 

reasonable value. For plasma contained within a radious r 

the number “ N of particles per unit length is given using (1) 

by N= r^n = lO^ n T*|^ / B and the power W reaching unit 

area of a wall of radius R is ( using (2)) 

3NkT 3X10^^ B T ^ -I 

” = TeT = 32 R-A ~ ' ^ 

8 

Again for a temperature of T= 10 iC/We find that the wall 
radius must be at least 

R=: 5.9 rv B^/^ (3) 

Now R must be larger than r and this condition 

yields/ 

B^3 ^auss (4) 
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For the minimum value of B when r is determined 
from above. 

A reactor with R=100cm is already a very large 

2 

device. Adopting this value and B = 2x 10 kG xvhich is 
practicable, we see from (3) that must be at least 1800 to 
satisfy criteria laid in (3) . With this value and E =0,01 
B must be larger than 180 kG which is consistent with our 
assumption of 200kG. Such values of have hot been achieved 
so far, even at temperatures smaller than thermonuclear 
onees but it may well be possible to reach this goal. We 
have presented these estimates to show that the search 
for alternative methods for obtaining fusion power is still 
necessary. 


The foregoing data imply confinement times of 
8 3 

Ims and densities of 10 /cm not easy to maintain for 
sufficiently long times for extracting power. 

This has led to an alternative scheme of confinement 
namely inertial confinement by lasers. It is for this reason^ 
that microexplosions involving much higher densities during 


much smaller times are tried using lasers. We now terja to the 
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problem of heating and compressing a pellet containing 
a deuterium-tritum misture. We see it now as a spherical 
object with a radios of the order of 100 micrometers . The 
scheme is to illimunate this target as uniformcly as possible 
xvith intense laser light whereupon some light will be 
absorlacd and materials will evaporate from the pellet 
surface. At high temperatures it will ionize and form a 
rapidly expanding plasma mantle. Under intense illumination, 
the plasma density may rise to such an extent as t o make it 
opaque to laser light. The heat from the absorption region 
will be conducted to the pellet where it will cause 
further evaporation. The pressure of the plasma fio ' •iUfc'med 
may become very high that a shock wave is launches into the 
target. The mass loss of the pellet by fast evaporation is 
called ablation. 

Thus according to our criteria, with E = 0.01, 

It 

nT= 10 , an explosion lasting lOOps would correspond 

to a density of 10 /cm which would be just sufficient 
provided that the fuel could be kept at this density 


during this time 
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Following is a brief review of the various 
fusion devices being experimented on. 

THE PINCH Device 

The theory of the constructed gas current was 
first developed by W.H.Bennet in 1934. A few years later 
a different approach was presented independently by L.Tonks 
who used the term "pinch effect" in the sense in which 
it is employed at the present time. Subsequently , a number of 
other theoretical analyses have been published of the 
construction of a current filament due to the action of the 
azlmaichal magnetic field generated by the current itself. 

Using the Benet relation I^= 2N k (T +T. )c^ it appears that 
the attainment of a sufficiently high temperature in a pinch 
discharge depends on the ablity to pass current of that 
magnitude. 

Since the temperature in the quasi-equalibr ium pinch 
varies as the square of the current*. IBalt equilibriufli' 
temperature is limited as the rate of energy loss will increase 
rapidly with the pinched discharge current. 
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A considerable amount of work has been done with 
"fast discharges" in which the constricted plasma is 
established in a time that is short in comparasion v/ith 
the sound transit and the magnetic field diffusion timeg. 

The dynamic pinch ( or shock pinch) as described above 
has received a great deal of attention experimentally since 
it offers^, in principle, a way to obtain high average particle 
energy in a plasma. When the pinch reaches its smallest radius 
the kinetic energy of the particles in the plasma reach their 
maximum inward velocity. If the pinched discharge column 
can be kept contracted long enough to permit randomisation,, 
i.e., thermal izat ion of the ion velocity to take place ,a 
highc. temperature will result. 

One of the major probems associated with the 
constricted discharge and infact, with any system of 
charged particles confined by a magnetic field- is tbat of 
instability. There are two main types of instablity which 
leads to the destruction of a pinched plasma. The first 
often called the sausage type of instability, results from an 

azimuthally symmetrical localized constrjjiition (and expansion) 
of the plasma column. 
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This is diagramatically represented in figure!. 1 
Consequently, the constrictions tend to become more and more 
constrictea , whereas the bulges will grow larger and larger 
until the pinched discharge is completely disrupted. 

The second type called toink instability originates 
with toe formation of a band or a kink in the plasma column 
while the latter retains its uniform ^ circular cross-sefction. 

This is diagramatically represented in figure 
It is seen that tne lines of force of the azimuthal self- 
magnet icfield , due to the current in the plasma, are brought 
closer together on the inside, but they are further ap^rt on 
the out side of a bend. Thus once a slight kink developes ,it 
will grow in size until the pinched dischaipge strikes the walls 
and is cooled. The 4^r plasma will then become diffuse, and fills 
the containing tube. 

A part from the problem of instabilities, the theta 
pinch devices are also bogged by problems of end losses. This 
has le-d to the the development of toroidal devices. 
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The adiabatic in variance of th-- magnetic moment in a 
space dependent mganetic field leads to a magnetic field ‘ 

configuration kno’im as a magnetic mirrors. * 

Tne other concept led, to development of a toroidal 
pinch device initially. The most promising of all thermonuclear 
reactors to date is of the Tokamak series developed in the 
Soviet Union, This is a toroidal pinch v;ith a very strong 
toroidal magnetic field generated by the current A poloidal 

field results from a large current induced in the plasma by an 
iron core. The stellarator was initially conceived as a toroidal 
like steady state confinement. The simplest form is that of a 
system of coils xvhich generates a magnetic field in a toroidal 
region of space. By contradistinction with the toroidal pinch 
the magnetic line of force here are "parallels" of the 
tor a id a .instead of mediams. In addition there are a number of 
other devices like the spheromak and base-feall machine. 

Recently ^ ir. has been found that in order to achieve a 
practical fusion device one has to use a high beta fusion 
device. Given the low^beta of tokamaks, there is renewed 
interest in theta pinch devices and more recently the 


REVERED =.PIELD PINCH (RPP) 
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The main plasma paramett.r of a thermonuclear 
configuration are determined by physical and economical 
consideration. For the output in Therma&nuclear energy from 
fusion reactor to exceed the energy needed to heat the plasma 
the -Droduct of the plasma density and the time T should 
satisfy the Lawson ceritcria; nT = 10 2 /cm v;here T is 

the mean time a charged particle can be held in the plasma 
volume. The mean pov^/er density desired from a fusion reacter 
is determined by the condition that the capital sost of the 
reactor should be comparable with that of the conventional 
power sources. The maximum power density possible is limited 
by the capability of the srounding material to withstand the 
energy released, xvithout being destroyed. Under these two 
condition, the power density released by the fusion reacter 

3 

should be btween 1 and 10CS''!'/cm . It. is clear that power 
release will be a function of the temperature and the density 
of the confined plasma. Under the circumstances, it is possible 
to consider three regions of confinement: these are: 

dO high density low temperature, short life tine 
plasmas. 
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(2) lovi density, low temperature, long life time 
plasmas . 

(3) low density, high temperature, lone life time 
plasmas . 

Whatever is the form of the magnetic confinement 
che plasma can always escape from the configuration. Energy 
losses from a controlled -thermonuclear reaction systems may 
be d^>vided into two broad categories- direct losses and 
indirect losses. Among the direct losses are those due to 
plasma instabilities, radiation, charge exenange, thermal 
conduction and diffusion. These are all energy losses that 
take place directly from the plasma itself- The indirect 
losses are those arising from esxternal circumstances, such 
as the losses resulting from joule heating in the magnetic 
field coils and from the conversion of the thermonuclear 
reaction energy into useful power. 

Some of the parameters of a thermonuclear reacting 
system, notably the dimensions of the reacting chamber and the 
stre^ngth of the confining magnetic field, can be yar ied 
in such a manner as to increase the power ^ output and 
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minimise losses. The broad rules which are ^plicable are 
refered to as scaling laws. . 

One of the main uses of a zero dimensional 
simulation program such as this might be to predict scaling 
laws for magnetic confinement devices, 

1.2 SUMS JSE JJL DI:FINT2;0NS ; 

In this section/ we first give definitions of system ■ 
parameters like beta, neutral pressure , crowbarring, elastic 
and inelastic collusion, ion an mean free path, and ambipolar 
diffusion. 

BETA:- It is often conveitjtient to determine the 
ratio of the kinetic pressure of the particle to the 
external magnetic pressure. The dimensionless quantity 
is then defined by 

' BV8/T 
o 


The energy balance equation is given by 


P-i- 

stt 
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This may be written as 
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Since the minimum value of B is zero, the ratio 
has a maximum value of unity; this I’Jould represent the 
ideal case of a perfectly diamagnetic plasma from which 
the magnetic field was completely occluded. In this event 


we have 


P, 


- /B 


max 


2 
o 

8;r 


where is the maximum kinetic pressure of a plasma that 

can be co'fined, in a study state, by a external magnetic 
K 

field of strength 

NEUTRi^L PRESSURE This is the pressure exerted 
by -che neutral gasf in the system. The magnetic region in 
vjhich the injected ion are trapped contain a large number of 
neutral gas atoms, owing to the finite initial pressure 
the system must have. 

CROWBARRIN3 This is a technique# by which 

the external magnetic field applied to the plasma is 

forcefully brought to zero after a predetermined amount of time. 

ION-ION MEAN FREE PATH:- This is the mean distance 

travelled by the ion betwec n c duQks ions with oter ion. 
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PLASI'IA ^CONFINEHENT TIME : This is the average 
time the plasma stays confined before it escape from the 
magnetic field. This is also equal to 1/e times the plasma 
decay time, 

AMBIPOLaR DIFFUSI0>J: This is a diffusion process 
which takes place due to the movement of electron and ions 
across the magnetic field with equal velocities. 

1.3 LITERATURE REVIEW: 

Computer simulation of various fusion devices 
using time dependent zero dimensional computer models 
have been carried out by several authors 

In reference (d^) the author has indicated that 
the modC^that he has developed was afest , simple^ time 
dependent , zero dimensional computer model for magnetically 
confine# plasmas. According to^author/. this method could used 
to model other devices or heating methods very easily. 
According to the author most of the terms do not depend 
on the particular geometry, containing at most the total 
plasma volume , and so could be used unchanged in modelling 


other devices 
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Several theoretical studies which describe and 
losses from theta pinch plasma column have been reported 
In the zero dimensional code of Green et al ( ) 

the electron and ion temperatures were assumed equal and 
particle end losses was included resulting in a solution 
for temperature which was seperated in spatial&time varialbles. 
The zero dimensional mod€4^ developed by Kievans et al 
includes both uniform plasma profiles as v/ell as the effects 
of ion thermal conduction and finite thermal conductivity. 

In the present thesis work, attempt has been made 
to include the geometrical th€ta pinch device as indicated 
in reference (2- ) into the model developed in refer eQce( ± ) 
to produce a zero-dimensional model of a pinch device 
where the preionzation has been done using a plasma gun, 

1.4 OUT LINE OF THE PRESENT WORK 

In chapter 2 a zero dimensional model is developed by 
defining the various energy loss terms. These energy loss 
terms have been space averaged in order to make them 
space independent. From this a set of equations, one each 
for plasma density electron energy density and ion energy 
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density have bc-en defined* In addition in chapter 2 the 
energy terms of a theta pinch device have been derived. 

In chapter 3 we established the method by which set of 
equations are solved. We also determine .in chapter 4 by 
parameter adjustment^ the optimal neutral pressure for 
obtaining the appropriate variation of temperature. As will 
be seen latter ^ any one of the follov;ing three situations 
can occur" in a pinch device: 

(I) Progressively decreasing electron temperature 
variation whereas a ^ progressively increasing densiti^. 

(II) Progressively incereasing electron temperature 
variation whereas a progressively decreasing densities. 

(III) Gradually approached peaks of electron ion 
temperaurs and densities. 

These cases are obtained by proper adjustment of 
neutral pressures. Next the electron and ion temperature 
variation in case of adiabatic compression studied. 

Finally, the surrmary and the main conclusion of this 
work arc presented in chapter 5. 
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Ca^vPTER 2 


The computer program (Sli’dULT) calculates the time dependent^ 
spitially averaged density/ electron temperature and ion 
temperature in a cylindrical multipole device [X ]• 
program is zero- dimensional in the sense that it predicts only 
spatially averaged quantities. It has b^^en used mostly to study 
continuous wave microwave heated plasmas, but pulsed microwave 
and gun-in jocted plasmas can also be studied. 

The numerical method consists of solving three coupled, time 
dependent, non-linear, first order, or'^inary differential 
equations, one for the p article density, one for the electron 
energy density and one for the ion energy density. 


dn _ 3 n §_ri 3 n 

dt S t ionisation 3t diffusion 3 t obstacles 

(a) (b) (c) 

- (l 

3t field decay. 

(D) 

dU 3U 3U 3U 

e _ e - e _ e 

dt 3t microwaves 3t excitation 3 1 ion collision 


(e) (f; (g) 

3U 3 U 

e _ e 

3t bremsstrahlung ^ 3 1 synchroton 
(h) (I) 



thermal conduction 


( J) 


( 2 ) 


dU, 


dt 


3U. 3U. 

1 i 

3t electron collisions 3t charge exchange 
(K) ■ (L) 

_ 

3 t thermal conduction 

(m) 


( 3 ) 



The first equation (l) describes the electron and 
ion densities in terms of ionization^ diffusion/ obstacle 
losses and field decay. The second equation (2) describes 
the electron energy density in terms of microwave heating/ 
excitation energy/ electron— ion collision/ bremsstrahluil^ losses; 
synchrotron losses and thermal conduction losses. The third 
equation (3 ) describes the ion energy density in terms 
of ion-electron collisions/ charge exchange and thermal 
conduction losses* 
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The computer results could not be compared with experimental 
values because all the required parameters in the pinch device 
V7ere not available* 

The tliree equations are solved by the simplest method of 
successive iteration using the current value of each parameter. 
The electrons and ion-tanperatures are increnented according to 

AT = (AU - T An)/n 
a e e ' 

ATj_ = (AU^-Tj_An)/n 
An is change in n in At. 

The iteration step At can be adjusted to insure any desired 
degree of accuracy. The initial conditions (n,T^, T^^) can be 
set arbitrarily. 

In this study a time dependent magnetic field as specif'-ed 
below. 

B = B sinut O) = 200 /sec t = 10”^sec. 

o 

In this study microwave power pulse'^specified in [ Jl ] 
was us ed . 

The neutral hydrogen gas is assumed to consist of thermal 

(T = T ) molecules. The initial value of the neutral density 
w 

is 011 ©xpstrim^n tsl p 3 .ir 3 in©t©ir whicTn must Id© spGcifisd.* 
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Reference [j,] indicates that the neutral (hydrogen) density 
is calculated at each time step from the pressure (p) as follov^s* 


n = 322p exp[-3.07xl0 


-4 na -17/T 


12.88p exp[ -1.7x10"^ — 


The first term represents the thermal (0,025eV) neutrals and 
np^^ represents Franck Condon (7eV) neutrals. The exponential 
factors account for the finite penetration depth of neutrals 
into the plasma. The thermal neutrals are assumed to be at 
wall-toiaperature (T^^) . The e^cperinient is assumed to take place 
on a time scale short compared vrLth the pumping time so that 
no neutrals are lost to t^e vacuum pumps- VThen the neutral density 
is increasing in time, the rate of increase of each component 
is limited to 


^^TH 1 . 9x10“^ 

"dF~ ^ i ^TH 


^FC 

dt 



representing the transit time of a neutral from the edge to the 
center of the plasma. 

(a) lONIZATION-The ionization rate for a Maxwellian electron 
distribution in a cold neutral gas, can be written in a form 
given by Drawn [ref(^ )]• For hydrogen the electron neutral 
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particles migrate across B by a random vralk process. VThen an 
electron collides with a neutral atom, the electron leaves the 
collision in a differoit direction but its phase is changed 
discontinously. 


D 


In such a case, D . = 2 ~~ 2 " 

1+w^T 


When 1, 


KT 1 

■X 




.,2 2 
“c^ 


KTy 


moJ 


0) 


eB 

m" 


=> D = 


KTxn^xcTv X m 
o 

2 

mxe xB 


KTxn^ cv jon 

“72 2 
B^xe , 


T . = X n 

-KTxn^x^vxn ^ ^ 


Kxnn^T 
o e 

n2 2 

B xa 


Consider the I term • 

n xnxEkT vnT 

In a fully ionized plasma, ^ =-^[T is in eV] 

B B 

T) xn n 

Electron current t = 5 x T 

B 


Kxn. xnx n 



Since 



= constant times 




(]ct)-V2 
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collisional ionization rate is 

+ ln(l.25(ltl/S) )) . 
where S = 15,6/T . 

^ Cl 


dn 

dt ionization 


- 3 71e~^ 
1+S T 

e 


1 _ 

20+S 


From another approach the ionization rate for a Maxwellian 
electron distribution is approximated by an anlytical fit 
for the data given in [n]. In viev/ of the inavailibility of 
ref[ ii. ] a analytic fit of the cross-section for ref [ '7-] was 
performed and a similar variation was found. In general the 
ionization is largely proportional to the electron temperatures. 


, nn 

dn _ CD . P -IVt 

dt ionization” *^^1/8 ® 

e 

(b) DIPFUSIOr^- Classical diffusion due to electron-ion and 
electron neutral collisions is given by 


2 3 nn T 

^ Q n Q- 1 ^ 

dt diffusion ” *'^■^^2 2^1/2 „2 2 

B a T ' B a 

e 

(I) (II) 

The above term has been satisfactorily derived as follows. 

Consider the II term. In a magnetic field charged particles 
will move along B as 

= nV. = +ia.nE-DX„- 


- J 


J n 


If there were no collisions, particles would not diffuse 
at all in the perpendicular direction- When there are collisions 
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-3 /2 

constant x nxn^xT ' xT 

~ diff 12 r~ 

3 xa 

constant xnxn 
_ o 

B^xa^xT^^^ 

Here ambi polar electric field has been considered » 

(c) OBSTACLES I Hoop supports and probes which intercept 
the plasma are assumed to collect particles at a rate given 
by the ion saturation current for a isotropic plasma* 

Current density J = nx velocity 

= constantx ^(T^+T^) 

2 

A^/a is a factor which determines the fraction of the 
area occupied by the obstacles. 

an 2,0x10^ C-Ter'TtS 

Hence we have obstacles = E ^ 

3. Xj 

(b) FIEI-iD DECAY I When is negative field lines are 

leaving the machine, and plasma is assumed to leave also at a 
rate given [ 1 ] 


3 n 

a t 


field decay 


-min[0. 


(e) microwave HEATING I Ref[ \ ] 



calculate the microwave 


heating rate as 
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_ej 

3 t ' microwaves 


9 

2x10 R 
Cl+n ^n) a^L 
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where is the density above which the microwaves are 

totally absorbed. In general, according to ref [ \ ] n^ is 

2 

given very roughly by n^ = 487f /Qq- 

Here is a brief discussion about the same. The incidence 
of a microwave on a plasma gives rise to the following phenomena 
(l) Absorption (2^ Diffusion, reflection and secondary emission 

( 3 ) Phase lag, caused by differences in refractive indices. 

( 4 ) Modification of the polarization and birefringence among 
other phenomena. 

Here Q. factor is the measure of the sharpness of resonance. 

Instead of a quantitative derivation of the expression on 
the R.H.S., a dimensional analysis in order to understand the 
derivation of the various terms is given below. 

Consider R.H.S. We have ^ L.H.S. 

cm'^ X ( f actor ) tim exarea 


A little more meaning could be given if one considers the 
process to be due to transit time magnetic pumping. If it is 
assumed that there is equipartition of energy between ions and 
electrons, P = 3nh ^ ,‘ This assumes that the phase velocity 
G)/k for optimal heating is of the order of the mean ionthermal 
velocity. 


0.^2 


P = 


K 


dU 

dt 


Here K = 
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) EXCITATION ; An electxon can also exoerience a large 
number of different types of inelastic collinsions with atom, 
ions and molecules. In a collision with one of these heavy 
particles, an electron may merely excite the heavy particle. 

The cross-section will be dependent on the energy of relative 
motion and the particular quantum level excited. Rsference['X ] 
has assumed that there is a 30eV/collision.rIence we have 


, 

§ ^ _ 30 ^ 

3 t excitation dt ionization 

(G) ion collisions . Consider . encounter between two 
charged particles. 'ere each particle moves in a hyperbola 
relative to the m of the two particles. The distance of 
closest approach is p ''the impact rsrameter' ' . Consider a 
test particle with mass m, a charge Eq/c and a velocity W . 
Then AW represents the charge of velocity of the test particle- 


Consider the Maxwell -Boltzmann distribution 

f ^ t3 t 2,,2 

frr\ nL _-L W 
f UJ; = -g - jg" - e 

where n and m are the particle density and the mass of the 
particle in question, and L is defined as L^ = . 

To measure the rate of diffusion in the W^ direction, 
we set N equal to the average number of encounters in one 
second. The resultant value of (AW^) , maasxiring the 

increase of ./elocity dispersion of a group of particles/second 
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2 

may then be denoted by < (AW^) ^ . This quantity is known as 
* * diffusion coefficient' ' . 


When the distribution function the field particles 
is assumed to obey a Maxwell -Boltzmann distribution only 
three independent diffusions need to be considered, < ^ • 

<(AW^^) > , and < (AW )j, > . The first quantity represents 
the rate at which moving test particles are slowed down by 
interaction with the field particles. The quantity < ^'^^11^^^ 
represent the rate of increase of (AV7) in a direction parallel 
to the origional motion of the particles. The corresponding 
quantity < (AW )j_ > represents the rate of increase in 
the perpendicular direction 


= h/p^ where h is the Debye Shielding distance and 
Pq the impact parameter 

3 /k^T^ _^l/2 


2ZZ, ^e 


As In increases T increases,* InAiocreases as n increases. 


Detailed computation of the diffusion coefficients have beai 
carried out by Chandrasekhar [ 3 ]. The resultant formulas 
for the three diffusion coefficients are 

<AW3_i> = ^)G(L,w} 

<(AW^j_)^>= ^G(L,W) 

A 

<(AW )^> = ~ (L,w)-G(L,w) 

W 
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where the '^diffusion constant'' Aj^ is defined by 


8 /I e^n^^Z^zJln 


m 


^(x) is the usual error function. 


(f> (x) 


^ 1/2 
71 ' O 


X 2 

/ e-y dy 


and the function G(x) is defined in terms of ^ (x) by the 
relationship 


G (x) = 


<P (x)~x (x) 

O 2 

2x 


The ''time of relaxation" denotes the time in which 
collisions produce a large alteration in some original velocity 
distribution. If one considers certain times which are defined 
in terms of the diffusion coefficients, we have ' t * 
is the time between collision in which reflection gradually 
deflect the test particles by 90°. In other words 

< (AW )^ > tj^ = W^. 

An energy exchange term t^ may also be defined by the relation 

< (Ae)^> tg = 

The change of energy, AE in a single encounter is given by 

AE = § { 2WAW^+(AW^^)^+(AW )^} 

Spitzer has determined the rate at ^irhich equipartition of 
energy is established between two groups of particles having 
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The actual calculation is difficult because V and V are not 
vrell knovm. Hence we assume that on electron of speed v has 
an impact parameter b with respect to an ion of charge 
which is assumed to be infinitely massive. The electron 
acceleration is 


. 

V ^ 


4 71 t mb 


2 


with V << C and neglecting the cross-product, we have 


dU 

dt 


2.2 
q V 
^e 


6 n £ C~ 

o 


dt 


2 4 
q. a 
-e 


96Tr^£^C^m^b‘^ 

O 3 


The total energy loss rata by the electron colliding with ions 
at all impact pararaeters is 


dt 


q^q'^Cn. V ) 
re 

96 

o e 


max 


/ I - 


b . b"^ 

min 


Taking b 


max 


~ ^min = 


h 


2 7T mV 
e e 


W = 
X 


2 2 

qT q n . n 

24 n eic^m h 
^ o 0 


^ 8kT 


m 


48x10”^'^ Z^n.n 


Wn t'-i- /rm'' 
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Maxwellian velocity distributions, but different kinetic 
temperature. T and Next he used (a) and averages this 

over a Maxwellian velocity distribution. VJe have 


dT 

dt 


T^-T 

■^eq 


where t is the time of equipartition given by 


3mm^k 


3/2 


8(2n)^/^n, e'^ln ^ ^ 


( 1 + ^)V2 

1 


This implies 


Now 


(T^-T^)x a(27l)^/^n^Z^Z^e^ln 

" 3mm vV2( I + !2)3/2 


m m., 


dt “ dt 


dU 

j 

dt 


Kxn x(T -T. ) 
e i 


T e.) 


( ~ + 
m m 


3/2 


dU 


which compares well with 


61n (T -T. ) 
e 1 

^2/3 


(h) BR EM S S TR AHLUNG ^ According to classical electromagnetic 
theory, the total power U radiated from a charge q moving- 
with a velocity V and acceleration v is 


^ v^-(VxV)^/C ^ 

671 e (1-v^/C^)^ 

o 

.V 

*»* 0 % 


8781 ? 




SYNCHROIRO^f RADIATION T Again according to basic 


electromagnetic theory. 


dt 


6 71 eC'^ 
o 


2 

In this case (v/c) << 1 


-(Vxv)^/C^ 

and V = . 


The electrons are slightly relativistic* Their angular 
frequency of rotation is determined by their relativistic 
mass, so that 






m 


I , I (■ 1 2 ,~2 \ 1/2 

=1 (.1-v /C ; ' 


Substiuting v^^ 


WqV we g et 


dU _ 
dt 


qV b2 


6 71 /C^ ) 


1 

J 


Expressed in terms of momentum. 


dt “ 


q^B^ 

e 


6 TC e m'^C^ 
o 


t 


which is to be integrated over the relativistic distribution 
f(p). Since f(p) is isotropic and since these are two direction 
associated with p but only one associated with p^^i. * ’”'® 
have 


2 


-2 
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The exact answer appears in terms of the raodified Bessel's 
function (z) of the third order. The integral definition 
of the Han<^l functions (z) of first kind and second 

order and the derviative relations for these functions 
are used to evaluate these integrals 

4„2 n kT K_ (mc^/kT ) 

„ = , _3 ■ a_ 

3 TC me K^(mc /kT ) 


From the asymptotic expansion of the modified Bessel's 
function^ 


^'^c "" 3" 7T £^ c 


n kT _ kT 

( c e ^ . 5 


-) (1+ I 


= 6.2xlO"^^B^n T [1+(T /204 ) 1 I'^atts ^cm^ 
e e 

Here is expressed in terms of B. 

(j) THERMAL CONDUCTION \ Heat conduction due to classical 
collisions and to obstacle losses [ 7 ] is given by 


dU ^ „ 

— -SS _ 2.5T ^1 

3 t thermal conduction e 3 t diffusion 


4- 


^.54T 


3 n 

e3 t obstacle 


(k) ^t.eCTRON collisions I Ions are heated by classical 
collisions with electrons. 


dU, 


dU 


dt electron collisions dt ion collisions 
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E EXCHAI^TG E I Reference [1 ] determined an 
approximate analytic form for change exchange losses, obtained 
by averaging cross-section found in ref. [ i2-] 

3U. 0 . 0186 nn 

1 _ o 1 

3 1 “ T ^+100 

Our first study had been to determine the exact expression 
itself. Later a sixth order polynomial was fitted for the 
data in ref [ />] and found similar variation in the cross 
section . 


(m) ion thermal conduction : 

collisions and to obstacle losses 


Heat conduction due to 
[ 7 ] is given by 


3U. 

i. I 

8 t thermal conduction 


2 . 5 T, 


8 n 
,Tt 


! + 2 T 1 

diffusion '^i 8 t 


obstacles 


first equation (1) describes the ■^l^ctron and ion 
densities ih'-terms of ionization, diffusion, obstapJt^ losses 
and field decay ."'"'^'l^ie-vSecond equation ( 2 ) ^jesrcr^^^s the 
electron energy density in'be^s o^-^rrticrowave heating, 
excitation energy, olectron;5.i6n cdl^l4sion, bremsstrahlung 
losses, synchrotron Ip^ses and thermal cohdi4ction losses. 

The third aquaj 

in terms ion- electron collisions, charge exchange aiM>thermal 

conddctinn losses. 


The above formulation of ref. [ 1 .] was modified in 2.2 
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to include plasma contraction terms as proposed by ref.[2-']* 
The following terms (n) and (P) contribute to the electron 
energy density term, whereas the term (0/ contributes to 
the ion energy density term. 

2 • 2 ADDITIONAL TERI'4S FOR A THETA - FINCH 1 

(n) a certain amount of work is done on the plasma by 

the system when the plasma contracts contributing to the 

increase in electron and ion temperatures. The following 

dA 

derivation is done for the term P %r = nT . 

s dt s dt 

For the theta- pinch device described in Chapter 1 we 
have the radial pressure balance as 


B^(r) = B?(r)+ 4 .03xl0“^°n(r) (T +T. ) 
e 1 e 1 

Here we have assumed that the radial inesrtial effects 

can be neglected and that the magnetic field has a small 

B 

curvature. Let 

a g- 
e 

Assuming a uniform magnetic field profile with 


B 


B(r) = 


r < a 


B / r > a 
e 


, ^ f 2.48 x10^^[b2/(t +T. )](H-S?) r< a 
n(ri=\ ‘-e-ei-' i ™ rh 

^ 0 r> a 

According to ref. [ ^] n = n(r) for a uniform profile. 
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3^ 

n - 2.A8xlO^^(^^(l^ 2) . 

3<p ^ ^2 1 

Using tha expression = 2 n a ,i _ 

a |j- d r r=:a 

o 


where 


<P = f B. (r)dp. 


Combining aquation [ ] and [ thu desired expression 

di 


fo^ :rr~ is obtained 


dt 


d^. 

c 

dt 


= l.( 


, dA T d3- 

1-- 4- i 


a dt ' B dt 

p e 


From another approach it can be established that 

_ i-Sf n dB dA ^ Y 

— ia ’.rh^re cc— ^ ^ -I ^ — T T (— i— — ) 

dt 2S -m^re g 3 A^ dt T +T.'- e 2 

a o p e 1 Tj 


'th 


T '^1 

-I- T, )h(0.25L-A.. , )+ - “ — 


1 1 
Tj 


hi 3 Tp e i 


From the above tv/o equations, for a uniform profile, we 
have 


dA, 


1 




dB (1-S^) 


+ 


^p (5-t^ ®e ^o'^i 

3. 


-A) 


where A 


= J ™(£ t-£.)+U^(-^= — )+T. (-T— )h(0.25L-X. . ) 


3 T-. G 1 ' B 


'th 


1 1 
T 


iX 


th 


(O) A terms similar to (m) above is used in the ion energy 
density expression. 
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( P) OHMIC HEATING I 

The ohmic heating tarn = L / ri (r)dA. 


With the assumption that the plasma is fully ionized the 
plasma resistivity 


h = 3.27xlO“^ln /T^/^[Ref(l) ] 


i 

j (r) from Ampere's law is 

3 


:i^ 

/j-. 


3 B(r) 
~Tr 


r) B (it ) 

For a uniform distribution ^ 

a r 


3 r 


[B (0.5+0.5(~)^)] 
e a -J 


0 . 5 B x 2 r B xr 
_ e _ e 

” ■" ”~ 2 “ - 2 
a a 

^ r -.2 2 

P Q "a /o B xr 

= L / 3 . 27 xlO"^x 2 . 2 /Ty ^ x x 2 n x r x dr 
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CHAPTER -=-3 

PROGRAM DESCRIPTION 

3.1 The numerical method consists of solving three 
coupled, time dependent, non“linear, first order, ordinary 
differential equations, one for the particle density, one 
for the electron energy density and one for the ion energy 
density. 

The three equation are solved by the simplest method 
of succesive iteration using the current value of each parameter. 
The electron and ion temperatures are incaremented ac«ording 
toj: 

T = (^ U -T A n )/n 
e e e 

T^ = (a- a n )/n 

The iteration step t can be adjusted insure any 
desired accuracey. The inital condition (n,Tg,T^) 
are given as input and can be set arbitrarily. 

3.2 The calculation is done in the following manner. 
Intially we specify the peak microwave PO in wats , UMAX 
the number of iterations, DEA the initial neutral density. 



TABLE i 

LIST OF DEVICES SIMULATED 


BIGLEV BIGUP 


TMAX 

0.043 

0.043 

A ( cm^ ) 

50.0 

50.0 

AL 

800.0 

800.0 

AO 

0.0 

700.0 

QO 

20000.0 

20000.0 


LITTLE 

QUADX 

PINCH 

0.005 

0.002 

0.005 

18.0 

6.0 

18.0 

270.0 

160.0 

800 

90.0 

3.0 

90.0 

2000 

500 

2000 


39 


TEA the initial electron temperature ^ TIA the initial 
ion temperature and TMAX the duration of the experiment in 
seconds. We also specify the plasma radijjus , length of the 
plasma, the obstacle area, magnetic field b 6, microwave 
cavity Q and the microwave frequencey F in GHZ, We specify 
the time step DT, the maximum neutral density DEMAX, the 
maximum electron temperature TEMAX, the maximum ion temperature 
TIMAX, thci minimum value of neutral density DNM, the maximum 
value of the magnetic field the microwave power , the ion 
saturation current and the maximum value of ereergy loss. 

We calculate the magnetic field and its derivative 
using the appropriate formula. The energy loss , BETA , the 
change in the neutral density DNT, are also calculated. After 
this the above calculations the current values of time DEA, 
TIA, B, P, Siiri,ELOSS §: BETA, are stored. 

Next the densities of electrons, ions and neutral 
atoms are incremented using the first order method. Following 
this WQ check for the time taken and the program is allowed 
to run or terminate appropriatc.ly • 
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it to sustain a continous reaction. (3) The energy 
released must be tapped in a controled manner. Now 
precisely the same cord ition arc needed to make a 
nuclear fusion reaction go and do useful work. The great 
difference is that for a fusion reaction the ignition 
point is rather high- hundreds of milli&ns of degrees 
cent^-rade . 

This one condition “ the artainment of which 
was quite unthinkable on the earth until recently - 
underlies all our problems. The. answer to the question 
"what is a controlled fusion reactor" is best given today by 
merely saying that it is a device within which appropriate 
isotopes of light elements could be caused to undergo 
nuclear fus^ion/ the end result being the controlled 
production and extraction of useful quanitities of 
energy^ in excess of that required to operate the device. 

Among the nuclear reactions which appear promising 
for a controlled fusion reaction are those involving the 
various isotopes of hydrogen, helium and lithium. Typical 
reaction are given Table (1). But unfortunately nuclear 
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3.3 SYSTEM PARAMETER ADJUSTMENT 

MAGNEiIC FIELD:- Any time dependent magnetic 
field can be specified. In this study v/e use 

B(t) = B^e ^^sincbt 

MICROWAVE POWER:- Any tim.e dependent microwave 
power pulse can be specified. Reference ( i. ) has used, 

P = MAX (O, 2B/B^-1) 

Though the initial conditions could be set quite 
arbitrarily, . . in this study the initial conditions have 

been set according to references ( X- ) • A very interesting 
outcome of the simulation is regarding the independence of 
the plasma parameter variation with initial electron and 
ion temperatures on the one hand and their sensitivity the 
neutral pressures on the other. 

3.4 ADDITIONAL PROGRAM FEATURES:- An inter actiTze graphics 
package has been incorporated in the program to obviate the 
time consuming excercise of plotting the parameters on paper. 
In addition to the large amounts of time saved, it also 

ensures greater accuracy in f lotting. 

The program is adequately representated in the flow 
chart given in Figure 
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CHA.PTER~4 

4.1 This study simulates a number of multipole devices and a 
theta-pinch device whose dimensions are given in table [3-1]. 

The first four devices correspond to multipole machines 
simulated by the V-Tisconsin group [ i ] . The pinch devices 
is similar to SCYLAC [ 2 ^]. 

The intent of this chapter is to simulate the time 
variation of various plasma parameters like the temperature, 
density and Beta. With the help of the zero-di, ensional model 
of the magnetic confinement system obtained in Chapter 2, 
a detailed dynamic analysis of the system for various neutral 
pressures has been carried out. The main theme of this 
process is the prediction of the plasma scaling laws and 
determination of the parameter space in which the system 
may be made to operate. It is important here to recall the 
fundamental assumption associated with the zero-dimensional 
model, that all the variables in the computation can be space- 
averaged over the cross-section. The second assumption is that 
the resistivity of the plasma can be computed with the assumpt- 
ion that the plasma is fully ionized. This can be somewhat 
Justified as one finds that the density of neutral atoms 
remains at a constant initial value mhough the ion and electron 
densities grow to 10-15 times the neutral density- The third 
assumption is that the plasma confincament time Xp is a variable 
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dependent upon the interaction of neutral atoms with ions 
and electrons. The question of how correct are the above 
assumption is difficult to answer with any generality, and one 
is usually reduced to verifying the validity of the assumpt- 
ions in various specific.- The sacrific^in curred in the 

A. 

process thus limits its validity to a certain range of operating 
conditions. The most important asoect of these numerical 
studies is to determine whether all the experimental observations 
can be verified using the existing theory. Por this we analyse 
the zero dimensional model by numercally solving the set of 
linear equations for predicting the transion-t behaviour. 

This qualitative comparision helps us in confirming the various 
types o f processes which are assumed to take place. 

4.2 THETA- PINCH .‘ 

Pig4!4,2 ) is a plot of the time variation of the electron 
ion densities. The time variation of the electron or ion 
densities indicates that higher values of densities of electrons 
and ions are obtained at higher neutral pressures. This is due 
to large contributions of the ionization term. One also notices 
that the peak keeps on broadening, indicating a shift towards 
temperature equilibrium. There is also a noticeable shift 
in the peak or a delayed peaking with decreased neutral density 
indicating that the ionization term becomes effective only 
tov;ards the caid with the decrease of the loss terms. 
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Fig. 14,3^ Is the time variation of electron tanperatures. 

The time variation of electron temperatures indiciates the 
strikingly large values of temperatures that can be obtained 
using the pinch effect. The temperature variation also 
indicates a peak broadening indicating an early temperature 
equalization. This can be attributed the smaller number of 
neutral atoms present initially. The shifting of the occur ance 
of the peak towards the starting point indicates the predominance 
of electron energy loss terms with increasing neutral pressures- 

Fig* ( 4 . 4 ) Is the plot of time variation of ion temperatures. 
One observation which one can make on seeing the time variation 
of ion temperatures is that the was do not gain significant 
energy (_ 60 e”) . This is primarily due to the fact that the 
electrons being more mobile, absorb all the in coming microwave 
power and there is little time for ions to gain energy via elect- 
ron ion collisions. One also notices that the maximum ion 
toTiperatures occur at neutral pressures which are neither too 
high nor too low. Rather they occur around temperatures 
^ which are' optimum' as will be seen later on- 

Fig. ( 4 . 5 ) is the plot of the largest peak of 'BETA' in 
each of the simulation runs. Contrary to usual expectations, the 
BETA value goes tlirough a minimum somewhere around the 
'optimum' values of neutral pressures. The BETA values tend to 
be very large either for low neutral pressures or for high 
neutral pressures. 
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Fig. ( 4 . 6 ) is a plot the variation of the final plasma 
cross-section with various neutral pressure. The steep drop 
could be due to neumerical instability. 

Fig. ( 4 . 7 ) is a plot of the plasma cross-sectional area 
as a variation with the operating time. The gradual increase 
in the cross-section is due to increased electron and ion 
densities and temperature resulting in decreased mean free path. 
The decrease in the cross-sectional area is due to decreased 
electron energies towards the latt'cr half of the experiment. 

The final cross-sectional area is very law in sharp contrast 
to other values in the experiment prim.arily due to very law 
electron densities and tomp. 3 ratures . 

Upon gradually varying the values of neutral pressures 

12 

such that the initial neutral densities vary from 1.6 x 10 
11 3 

to 6 X 10 particles/cm it was found that one could obtain 
reasonable variations of densities and tanperatures only for 
neutral densities varying from 5.1x10^^ to 5~.8xlO^^ particles/cm^ . 
Upon assuming neutral densities either below 5.1x10^^ or above 

11 3 

5.8x10 particles/cm the density variation or the electron 

tanperature variation become insignificant. In this study 

11 3 

neutral densities of 531x10 particles/ cm was used as shown 
in fig. ( 4 . 2 )/ ( 4 , 3 )/ (4.4). 
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Fig. ( 4 . 8 ) is the plot of normalized values of plasma 
paramaters with time at optimum neutral densities (531x10^^). 

One notices an early and broad electron temperature peah. 

This is because the microwave power input is primarily absorbed 
by the electrons initially. As predictable the magnetic field 
variation is a sinusoid peaking somewhere midway. The electron 
or ion densities shows a peaking in the latter help of the 
sirnultation. The ion temperatures sho-t; a peaking midway 
between the electron and density oeaks . There is only a small 
space of time in which the ion temper ac ares, electron tenperat- 
ures and densities show a peak. Significantly though the elect- 
ron and density variations show a broad peak, the ion 
temperatures shows a sharp peak. The time lag of the B peak and 

density peak is noticable in the Wisconsin work ( d. ) . 

OCTO’Poie. 

Fig, ( 4 . 9 ^ is a plot of electron or ion density variation 
with time for various neutral densities. Unlike the theta- 
pinch , the density peaking occurs at the same location and 
there is apparently no shifting in the time of peaking. 

Secondly the peaking is also not as broad as that of the 
theta pinch. This is because the ionization term becomes 
predominant only midway in the simulation. 

Fig, (4,1Q1 is a plot of electron temperature variation 
with time for various neutral densities. One notices the 
increased maximum values of elactron temperatures with 
decreased neutral pressures. There is also a predominant 
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B 

Calculation of 2>^/dt 
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tail coming into the picture at low neutral densities. This 
giv ’S an indication that the losses which were predominantly 
due to Branmstrahlung and synchrotron were significant only 
then. There is also a considerable amount of peaTc broadening 
indicating that ideally the electron temperatures attains a psak 
value and retains it till the end of the simulation. 

Fig. (4.11) is a plot of ion temperatures varying 
with time. The ion gain temperature only in the latter half 
of the simulation indicating better aectron-ion coupling then. 
In fact in the last few tim.e stens the ion tomperaturas show 
an almost asymptotic increase in values. Even then the ions 
attain values of around (SOeV) at vrry high neutral densities. 

Fig. (4d2) is a plot of the variation of BETA with time 
at various initial neutral densities. It can be observed that 
maximum value of the BETA peak occurs at neutral densities 
which are around 17000x10^ particles/cm^ . This slight deviation 
is due to the fact that (nT) product maximum occurs at this 
neutral density. 

Similar to the theta“Pinch/ it was octermined that neutral 

12 1 

densities around 1.7x10 particles /cm were most suitable. 

Fig, (4.14) is a piot of normalized plasma parameters with 
time. Unlike the case of a theta pinch fig. (4.8) the electron 
temperatures and densities not have a broad peak* As one would 
have also noticed the maximum values of the various parameters 
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are several orders less than th@se corresponding^ the theta 
pinch* Also present in the ion and electron tanperature 
variation is the tail which attains very large values towards 
the end of the experiment. This is primarily due to reduced 
losses towards the end of the experiment. 

Ref. ( i' ) has indicated that the energy loss rate can be 
taken to be proportional to -dB/dt. Upon incorporating such 
term one can obtain the plot of energy loss rate versus time 
as shown in fig. (4J3) . Ref. ( i. ) has indiciated that this 
can be taken to be the rate at \vhich the magnetic field lines 
leave the system. Fig. (413) indicates that the energy loss 
rate grows rapidly with time reaching very large values towards 
the end of the experiment. There is no energy loss in the first 
half of the experiment because here the magnetic field lines 
grow stronger with time. 

THEORETICAL BASIS FOR NEUTRAL PRESSURE VARIATION t 

In this study the pinch formed after the pre-ionization. 
Although heating is not a problem for .the high energy injection 
device, a new difficulty arises. The magnetic region in 
which the injected ions are trapped contains a large number 
of neutral gas atoms, owing to -the finite initial pressure 
which any system must have. The neutral gas produces a steady 
loss of the trapped ions by the process of charge exchange. 

In this process the hot ion captures an electron from the 
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gas atom and thus escapes from the systgt as an energetic 
neutral atom. Thus charge exchange competes quite seriously 
with the desire'l nuclear reactions in the gas. 


It would thus seem that the first major barrier which must 
be overcome by a high energy injection device is the neutral 
gas background. 

DIRECT (BRUTE FORCE) PUI4PIN3 I 

The most obvious method, of course, is simply to reduce 
the initial pressure of the system until charge- exchange no 
longer competes with the usual plasma processes of Coulomb 
scattering and nuclear reactions. The loss rate of ions/unit 
volume can be shown to be limited by the charge exchange. 

NEUTR.giL ' BURHQUT" ^ 

It is possible to obtain a condition for creation of a 
thermonuclear plasma which is less restrictive than the one 
above. This takes advantage of the fact that the self 
shielding action of the plasma tends to lower the neutral 
atom density in the plasma interior. Neutral atcsns are removed 
in two ways. One is by the charge- exchange process itself 
since the fast neutral which leaves the plasma region. The 
second and more effective, process is simple ionization of 
the neutrals by the fast ions (and electxons) of the plasma* 

As long as the rate of destruction of neutrals in the 


TABLE 1 • 1 

TYPICAL FUSION REACTIONS 


1. D+D 

wm 

"He 

+n + 

3.25 

MeV 

2. D+D 

= 

T 

+p + 

4 

Mev 

3. D+T 

sc 

^He 

+n + 

17.6 

MeV 

4, ^He+D 

= 


+p + 

18.3 

MeV 

5 . ^Li+D 

= 

afne 

+ 

24.4 

MeV 

6 « "^Li+p 

== 

aiHe 

+ 

17.3 

MeV 
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plasma interior is small than the rate of influx of 
neutrals from external vacuum region/ there will be no 
appreciable difference between the neutral density internal to 
the plasma and external to the same. There is a critical 
point/ hov/ever for the injected current at which point the 
neutrals are being destroyed as soon as they enter the 
plasma region. Once this point is passed/ the neutral 
density drops in the interior and the plasma density increases, 
hence burningout more ' 'neutrals' ' . This critical point 
will bo called ''burnout'' . It is expected that the neutral 
density will drop sharply from its ^’■acuum value when the 
injected current is only slightly above the burnout point. 

A crude first estimate of the burnout current is 
obtainable directly from the considerations above. 

TRAPPED ION PUMPING (G^ITSRING) : 

A certain fraction of the neutral atomiS of the system 
become ionized upon passage through the trapped plasma region. 
These ions are then constrained to move along a field line 
and will drift out to the 'ends of the mirror container. If 
suitable pumping arrangements are then provided, this 
mechanism could represent an important addition to the total 
pumping speed. In some cases it may dominate and appreciably 
lower the entire pressure in the interior of the system. In 
this case the trapped ions, of course, 
as an ion- pump, or as a getter. 


are acting exactly 
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It has thus been demonstrated that several aspects of 
i_he dynamics of pinch formations can be brought out using 
even a zero dimensional model of the pinch device'. 
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CHAPTER-5 

5.1 SUMMARY AND CONCLUSIONS: 

In the present work we have developed a zeroed imensional 
model of a octopole and a chcta^pinch device for simulating the 
dynamic behaviour and properties of various plasma 
parameters. 

Chapter 1 is devoted to a general description and 
some useful definitions related to the present work. It also 
contains literature review and outline of the present work. 

In Chapter 2 we obtained zero-dimensional models of a 
theta-pinch device and a octopole by deriving the various 
energy loss or gain terms. Our aim was to study the time 
behaviour of various plasma parmaters under various initial 
operating conditions, such that one could determine if 
possible parameter spaces of feasible and stable operation. 

A Very acceptable variation plasma densities and 
temperature la is obtained at neutral densities between 
580 X 10® and 510 x 10® particles / cm® for a theta pinch. 

If the neutral dt-nsitles are higher than 580 x 10 ®partlcle!/cm® the 
variation of temperature s becomes very predominant with 
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the electron and ion densities tending to be very low. 

The variation of thi- plasma cross-sectional area is also 
very typical of a thcta-pinch device. The variation in 
temperature and densities for the octopole is several 
orders less compared to a theta-pinch device for a given 
set of initial conditions. 

5.2 FURTHER StCGPJSTFD WORK:- 

The first improvement that can be done cancerns 
the method adopted to solve the equations of the zero-r 
dimensional model. In this study a first order method has 
been used. More accurate results can be obtained if the 
equations were solved using a second order method (Appendix A). 
J.R, McGowan et al ( >^) have indicated in a paper that 
Heindrich et al found the rotational energy to 

be comparable to that of the thermal energy. Hence the 
code could be modified in order to be able to examine 
possible effects of rotation on transport phenomenanin 
the plasma. The zero D model includes five first order 
differential equations. Four of the equation are derived 
from expressions for the time rate of chance of electron and 
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ion tempi^ratures, total particle in ventory and the 
magnetic flux within the^. - plasma. The fifth equation is 
obtained from radial pressure balance. 

Another extension could be in the direction of 
developing a zero^-D model of a reversed fi^2dpinch. 
rier>_ in this type of fusion device/ the closed field lines 
arc produced internally, a progress report on "Study of 
plasma convection and wall interactions in magnetic 
confinement system" which isaUS department of energy 
project briefly explains the recently developed zero-D model. 
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APPENDIX A 


The three equations derived in Chapter 2 can be solved 
simultaneously by the following second order method* Given a 
time step dt and the time derivation h , T , T. at time t, then 
calculate n, at time t+dt. First calculate B, P^n^ at time 

t. Second increment n, T^, Tj_ to approximate them at t 40 .=; dt 
(x is any of n, 


x(t+0 »5dt) S5 x(t) + (0.5dt)(x ^ dt^ * 

Third# calculate at t+O.sdt# using values of at 

time t and n^/T^^Tj^ at t+O.Sdt. Fourth, calculate T^, T^ from 


i 


2 

3 


U 


e#i 


-T . n) /n. 
e, 1 


Finally, using t+0,5dt, increment the original 

n, Tg, Tj_ at time t to get (x is any of n, T^, T^) 
x(t+dt) = x(t) + Cdt) (x 


APPENDIX C 

Calculation of 

Using Maxwell's equations, 

V X B = J 

B 

V X £ = — 

Consider now the magnetic field diffusion eguation* The plasma 
magnetic field B^(r) is related to the external magnetic 
field, plasma column area, and resistivity by a magnetic field 



diffusion equation obtained frcm» Maxwell^ s equation and 
Ohm^ s law in the form 


Jn = ) 

where E is the azimuthal electric field and is the 
plasma radial velocity. Since the plasma diamagnetic current 
flows transverse to the magnetic field lines and since the 
magnetic field in the sheath is of sufficient magnitude 
so that the electron gyrofrequency is many orders of magnitude 
larger than the electron- ion collision frequency, the 
strong magnetic field limit expression for plasma resistivity 
is used in this analysis ] 

T) , =3.27 X 10”^ In * 

1 ^ e 

Reference [ ] indicates that the magnetic field equation 

can be expressed as 


= V X ( .t-L- V X S - U x 5) 

3 t 

With the long ccmpression oil length assumption, Br/Z <<1, 
the z-component can be written 


_ 1 ^ % 
~ r 3 r 3 ^ 


1 3 

r 3 r 


CrU^B^) 


Integrating the above equation over the plasma cross 

secUonal area and using the Leibnitz' s rule for differentiating 

integrals 
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Using A 
on both 


^ r and _ (g a/ 9t)(l/2H r), then the second terms 
sides cancel leaving 
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NOMENCLATUkE 


Vq = velocity of diffusion 
C = speed of light 
kT 5= tonperature in eV 
p = kinetic pressure 

= electron temperatures (eV) 

= ion-tenperature (eV) 
n = number of electron or ion-densities 
B = magnetic field in gauss 
= neutral density 
a =s plasma radius (cm) 

0 )^ = plasma frequency 

7?^ = resistivity 

L = length of plasma column 
= ion-ion mean free-path 

/ 9 3 

= average electron energy density VIO eV/cm 

7 3 

= average ion energy density ClO eV/cm'^ ) 

= obstacle area (cm ) 

P = microwave power (watts) 
f = microwave frequency (GHZ) 

Qq = unperturbed cavity of plasma chamber 

t = time (sec) 

s 

"^th" thermal conduction time of species 
Tp = plasma column particle confinement time 
Eg = average energy loss per species 
Tgq= electron energy equilibration term 



